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ABSTRACT 

Several  space-time  adaptive  processing  algorithms  have 
been  proposed  to  detect  weak  targets  in  the  presence  of  strong 
interference,  especially  clutter  and  jamming.  Except  for 
Displaced  Phase  Center  Array  (DPCA)  processing,  radar 
signal  processing  algorithms  ignore  the  fact  that  the  location 
of  the  clutter  ridge  in  angle-Doppler  space  is  known,  given 
the  platform  speed  and  direction.  This  paper  introduces  our 
attempt  to  exploit  this  a  priori  knowledge  in  conjunction  with 
the  joint  domain  localized  processing  algorithm.  Using  a  two- 
dimensional  filter,  clutter  is  rejected  in  a  first,  non-adaptive 
stage,  followed  by  adaptive  processing  in  the  angle-Doppler 
domain. 

1.  INTRODUCTION 

Airborne  surveillance  radar  systems  operate  in  a  severe 
and  dynamic  interference  environment.  Space-Time  Adaptive 
Processing  (STAP)  techniques  promise  to  be  the  best  means 
to  detect  weak  targets  in  such  interference.  STAP  algorithms 
process  data  in  two  dimensions,  separating  the  target  from 
interference  in  both  angle  and  Doppler.  Consider  a  phased 
array  antenna  with  N  spatial  channels,  and  M  pulses  per 
Coherent  Processing  Interval  (CPI).  The  most  straightforward 
STAP  algorithm  uses  all  NM  degrees  of  freedom  (DOF), 
estimating  the  NM  dimensional  covariance  matrix  of  the 
interference  to  minimize  the  expected  squared  error  with 
respect  to  the  desired  signal  [1].  It  is  now  well  accepted  that 
this  fully  adaptive  algorithm  is  impractical,  since  it  is 
impossible  to  obtain  enough  data  samples  to  estimate  such  a 
large  covariance  matrix,  and  the  associated  computation  load 
is  prohibitive. 

To  overcome  the  drawbacks  of  the  fully  adaptive 
algorithm,  researchers  have  proposed  several  algorithms  that 
limit  the  number  of  adaptive  weights  [2].  In  particular,  Wang 
and  Cai  introduced  the  Joint  Domain  Localized  (JDL) 
algorithm,  a  popular  post-Doppler,  beam-space  approach  that 
adaptively  processes  radar  data  after  transformation  to  the 
angle-Doppler  domain  [3].  Processing  is  restricted  to  a 
Localized  Processing  Region  (LPR)  in  the  transform  domain, 
significantly  reducing  the  number  of  DOF,  while  retaining 
maximal  gain  against  thermal  noise.  The  reduced  DOF  leads 
to  corresponding  reductions  in  required  sample  support  and 


computation  load. 

In  airborne  radar,  one  crucial,  indeed  often  dominant, 
component  of  the  interference  is  the  ground  clutter.  Due  to 
the  motion  of  the  radar  platform,  the  clutter  occupies  a  clutter 
ridge  in  the  angle-Doppler  domain,  i.e.,  after  transforming  the 
space-time  data  to  the  angle-Doppler  domain,  the  clutter  is 
localized  along  a  line  [2] .  The  slope  of  this  line  is  determined 
by  the  speed  of  the  airborne  platform,  a  parameter  that  may 
reasonably  be  assumed  known.  The  location  of  the  clutter  in 
angle-Doppler  space  is  therefore  known  a-priori.  This 
information  has  been  used  in  non-adaptive  processing,  most 
commonly  by  DPCA  processing,  however  has  been  ignored 
in  the  development  of  STAP  algorithms. 

A  central  question  asked  in  this  paper  is  how  to  exploit 
this  a  priori  information.  This  work  develops  a  simple 
extension  of  the  JDL  algorithm  that  uses  knowledge  of  the 
location  of  the  clutter  ridge  to  pre-filter  the  clutter.  The 
resulting  space-time  signal  is  then  processed  to  detect  weak 
targets.  A  crucial  factor  in  the  success  of  this  algorithm  is  the 
ability  to  model  the  changes  in  the  space-time  steering  vector. 
The  model  for  the  change  in  the  space-time  steering  vector 
was  developed  in  [4].  With  minimal  increase  in  computation 
load,  the  pre-filtering  approach  results  in  significant 
performance  gains. 

This  paper  is  organized  as  follows.  Section  2  briefly 
reviews  the  JDL  algorithm  as  modified  in  [3]  and  then 
introduces  the  pre-filtering  technique  developed  here.  Section 
3  presents  numerical  examples,  using  both  simulated  and 
measured  data,  to  illustrate  the  performance  improvements  to 
be  had  by  using  pre-filtering.  Section  4  ends  the  paper  with  a 
summary  and  some  conclusions. 

2.  CLUTTER  PRE-FILTERING 

We  begin  with  a  brief  summary  of  the  JDL  algorithm  as 
developed  in  [4].  In  the  JDL  algorithm,  target  and 
interference  are  transformed  to  7/a  angle  and  rjd  Doppler  bins 
forming  an  7/a  x  r/d  LPR.  This  is  achieved  using  a  trans¬ 
formation  matrix  comprising  the  space-time  steering  vectors 
corresponding  to  the  angle-Doppler  bins  in  the  LPR. 
Representing  as  a(0)  the  length-A  spatial  steering  vector 
associated  with  angle  0  and  b(  fd  )  as  the  length -M  temporal 
steering  vector  corresponding  to  Doppler  frequency  fd,  the 
transformation  matrix  T  for  7/a  =  rjd  =  3,  is  given  by 
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where  the  LPR  encompasses  angles  0_i,  0O  and  0i  and 
Doppler  frequencies /_i,/o  and/i  and  s(/j,  0)  =  b(fd  )  ®  a(0) 
represents  the  length-Mf  space-time  steering  vector. 

The  transform  domain  data  and  transform  domain 
steering  vector  are  given  by 
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where  the  tilde  (  ~  )  represents  the  transform  domain,  x 
represents  the  length-TVM  space-time  data  vector  and  s  the 
length-TVM  target  steering  vector,  x  represents  the  length- 
77a/7d  vector  of  data  within  the  LPR  and  w  the  length- %% 
adaptive  weight  vector  in  the  transform  domain.  The  matrix 
R  represents  the  estimated  covariance  matrix  in  the 
transform  domain. 

In  this  work,  the  space-time  data  is  first  filtered  using  a 
two  dimensional  band-stop  finite  impulse  response  (FIR) 
filter.  The  stop  band  of  this  filter  is  chosen  to  match  the 
location  of  the  clutter  ridge,  i.e.,  the  filter  rejects  the  clutter 
and  allows  all  other  signals  through.  The  filter  is  designed 
using  frequency  sampling  [4].  The  design  process  requires 
identification  of  the  filter  response  at  a  finite  number  of 
points  in  angle-Doppler  space.  The  frequency  response  of  the 
resulting  FIR  filter  is  then  guaranteed  to  equal  the  desired 
response  at  these  points  in  angle-Doppler  space.  The  process 
and  results  in  a  matrix,  W,  of  filter  coefficients. 

Later  in  this  paper  we  present  results  using  simulated  data 
wherein  the  clutter  ridge  has  unit  slope  in  the  spatial- 
frequency  -  Doppler  frequency  (angle-Doppler)  domain.  The 
ideal  two-dimensional  (2D)  filter  is  shown  in  Fig.  1  while  the 
practical  filter  used  in  this  work  is  shown  in  Fig.  2.  Clearly, 
the  width  of  the  band- stop  filter  is  an  important  design 
parameter.  However,  analysis  of  the  impact  of  this  parameter 
is  beyond  the  scope  of  this  paper. 

The  space-time  data,  x,  and  space-time  steering  vector  are 
filtered  using  the  band-stop  filter  coefficients. 


Figure  1:  Ideal  Two  Dimensional  Bandstop  Filter 


Figure  2:  Bandstop  filter  used  in  pre-filtering  process 

The  JDL  adaptive  process  continues  using  Eqns.  (2)-(4)  with 
the  filtered  data  x  replacing  the  original  space-time  data  x 
and  the  filtered  steering  vector  s  replacing  the  original 
space-time  steering  vector  s.  Note  that  incorporating  the  filter 
is  made  possible  only  using  the  transformation  matrix 
introduced  in  [3]. 


X  =  W*X,  (6) 

S  =  W*S,  (7) 

where  X  and  S  represent  the  space-time  data  and  steering 
vector  arranged  as  M  x  N  matrices,  *  2D  convolution  and  the 
hat  (A)  above  a  variable  represents  filtered  data.  The  new 
transformation  matrix  is  then  given  by 


T  =  [s(/_, ,  0_, )  s(/_, ,  0O )  s(/_,  ,0,)s(/o,  0_, ) 

s(/o,0o)s(/„,0i)  s(/j ,0_, )  s(/j ,0O)  s (/L0,) 


3.  NUMERICAL  EXAMPLES 

This  section  presents  numerical  results  illustrating  the 
gains  due  of  the  pre-filtering  scheme  described  in  Section  2. 

3.1  Simulated  Data 

The  first  example  uses  N=  12,  M  =  14  and  a  filter  size  of 
21  x  21.  The  filter  is  designed  to  approximate  the  ideal  filter 
in  Figure  1 .  Figure  2  plots  the  frequency  response  of  the  filter 
used  in  this  example.  Table  1  lists  the  parameters  of  the 
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antenna  array  and  interference  used  in  this  example.  The  table 
also  lists  the  parameters  used  in  the  processing. 

Figure  3  plots  the  angle-Doppler  response  of  the  data 
within  the  target  range  bin.  The  clutter  ridge,  with  slope  1,  is 
clearly  seen.  The  target  amplitude  is  artificially  enhanced  to 
be  visible. 


Target  range  bin  spectrum  (before  filtering) 


Normalized  Doppler 


Figure  3  :  Angle-Doppler  clutter  ridge  before  filtering 

Figure  4  plots  the  angle-Doppler  response  of  the  post- 
filtered  data  within  the  target  range  bin.  The  clutter  ridge  is 
clearly  eliminated.  Note  that  due  to  the  type  of  filter  used,  the 
target  information  does  not  appear  corrupted.  However,  by 
accounting  for  the  change  in  the  steering  vector  using  Eqn. 
(7),  any  impact  of  the  filter  on  the  target  is  accounted  for. 
This  accounting  for  the  change  in  steering  vector  is  also  taken 
into  account  in  the  JDL  algorithm  in  Eqn.  (3). 


Target  range  bin  spectrum  (after  filtering) 


Normalized  Doppler 


Figure  4:  Angle-Doppler  clutter  ridge  after  filtering 


Table  1:  Parameters  for  example  using  simulated  data 


Number  of  Pulses 

14 

Number  of  Elements  in  Array 

12 

Element  Spacing 

0.5  X 

Pulse  Repetition  Frequency 

1024  Hz 

Transmit  Array  Weighting 

Uniform 

Number  of  Clutter  Patches 

181 

Target  Signal-to-Noise  Ratio 
(per  element-pulse) 

OdB 

Normalized  Target  Spatial  Freq 

0.1 

Target  Normalized  Doppler 

0.4 

Thermal  Noise  Power 

Unity 

Clutter-to-Noise  Ratio 

5  OdB 

Number  of  Doppler  bins  in 
LPR 

3 

Spacing  between  Doppler  Bins 
(normalized  frequency) 

0.05 

Number  of  Angle  Bins  in  LPR 

3 

Spacing  between  Angle  Bins 
(normalized  frequency) 

0.05 

Number  of  Filter  Weights  in 
Spatial  Dimension 

21 

Number  of  Filter  Weights  in 
Time  Dimension 

21 

Normalized  Filter  Notch  Width 
(in  Angle-Doppler  space) 

0.14 

15  20  25  30  35 

Range  bin 


Figure  5:  Impact  of  pre-filtering  with  simulated  data 

Figure  5  plots  the  MSMI  statistic  versus  range  bin  using 
JDL  with  and  without  pre-filtering.  Figures  1  and  2  illustrate 
the  ideal  and  practical  filters  used  to  suppress  the  clutter 
ridge.  A  weak  target  (OdB)  is  inserted  in  range  bin  26.  Note 
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the  improved  discrimination  between  the  target  and 
surrounding  range  bins,  an  improvement  of  approximately 
1.5dB.  It  is  important  to  emphasize  that  this  gain  is  achieved 
purely  due  to  pre-filtering  the  clutter.  The  elimination  of  the 
clutter  before  translating  the  data  to  angle-Doppler  space 
implies  a  lower  interference  level  in  angle-Doppler  space, 
allowing  for  improved  performance.  The  next  example 
illustrates  the  performance  improvements  in  measured  data. 

3.2  MCARM Data 

This  example  illustrates  the  impact  of  pre-filtering  the 
clutter  when  using  measured  data  from  the  Multi-Channel 
Airborne  Measurements  (MCARM)  program  [5].  The 
MCARM  database  includes  a  vast  collection  of  clutter  and 
signal  measurements  collected  by  an  airborne  radar  over 
multiple  flights  with  multiple  acquisitions  on  each  flight.  The 
acquisitions  used  in  this  example  uses  a  N=  22  element 
rectangular  array  arranged  in  a  2x11  grid.  Each  CPI 
comprises  128  pulses  (M=  128).  Also  provided  with  the  data 
is  a  set  of  measured  spatial  steering  vectors  for  some 
specified  azimuth  and  elevation  angles.  This  example  uses 
data  from  acquisition  575  on  flight  5,  acquisition  cycle  ‘d’. 

This  example  uses  an  injected  target  with  amplitude 
chosen  such  that  it  cannot  be  distinguished  without  adaptive 
processing.  The  target  is  injected  at  broadside  and  Doppler 
bin  -9  into  range  bin  290. 

The  use  of  a  2x11  array  requires  careful  filtering  of  the 
MCARM  data.  The  notion  of  the  clutter  covering  a  single 
ridge  in  angle-Doppler  space  assumes  a  linear  antenna  array. 
The  data,  therefore,  is  filtered  in  two  batches  using  11 
elements  in  each  batch.  Furthermore,  the  use  of  measured 
data  implies  the  standard  steering  vector  (of  phase  shifts)  is 
not  valid  and  the  measured  steering  vectors  must  be  used  in 
the  pre-filtering  process.  Figure  6  plots  the  clutter  ridge  in  the 
MCARM  data.  The  amplitude  of  the  injected  is  increased  by 
20dB  to  be  visible  in  this  non-adaptive  processing.  The 
clutter  ridge  is  clearly  visible. 


Target  range  bin  spectrum  (no  filtering) 


Normalized  Doppler 

Figure  6:  Angle-Doppler  clutter  ridge  before  filtering 
(MCARM  Data) 


Figure  7  plots  the  angle-Doppler  clutter  ridge  after  using 
the  band-stop  filter.  The  clutter  ridge  is  clearly  eliminated 
without  significantly  affecting  the  target.  The  filter  uses  21 
coefficients  in  both  the  spatial  and  temporal  dimensions. 


Target  range  bin  spectrum  (with  filtering) 


Normalized  Doppler 


Figure  7:  Angle-Doppler  clutter  ridge  after  filtering 
(MCARM  Data) 


MSMI  with  no  filtering 


Figure  8:  MSMI  Statistic  versus  Doppler  and  Range 
(No  Pre-Filtering) 

Figure  8  plots  the  modified  sample  matrix  inversion 
(MSMI)  statistic  versus  Doppler  bin  and  range  without  any 
pre-filtering.  While  the  target  is  visible,  false  alarms  due  to 
the  clutter  are  also  visible.  Figure  9  plots  the  MSMI  statistic 
when  using  the  pre-filtered  data.  As  is  clear,  the  false  alarm 
rate  due  to  the  clutter  has  been  significantly  reduced.  The 
injected  target  is  clearly  visible.  Note  that  the  other  targets 
visible  in  this  figure  were  identified  earlier  as  corresponding 
to  highways  with  vehicles  at  highway  speeds  [6]. 
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Figure  9:  MSMI  versus  Doppler  and  Range  using  pre-filtered 

data 
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4.  CONCLUSIONS 

This  paper  has  presented  a  pre-filtering  approach  to 
eliminating  the  clutter  before  adaptive  processing.  The 
algorithm  proposed  here  uses  the  fact  that  the  clutter  lies, 
mainly,  on  a  single  ridge  in  angle-Doppler  space.  Design  and 
use  of  a  two-dimensional  band-stop  filter  eliminates  all 
signals  lying  along  this  ridge.  Numerical  examples  using  both 
simulated  and  measured  data  illustrate  the  performance  gains 
available  by  pre-filtering  the  clutter. 

The  approach  suggested  here  is  relatively  simple  and 
does  not  require  significant  additional  processing.  In  this 
regard,  it  is  hoped  that  the  work  here  is  a  first  step  to  more 
sophisticated  algorithms  that  exploit  a  priori  knowledge  in 
the  adaptive  process. 
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